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Abstract 

We analyze the possibility of heavy-light tetraquark bound states by means of 
a chiral constituent quark model. The study is done in a variational approach. 
Special attention is paid to the contribution given by the different terms of the 
interacting potential and also to the role played by the different color channels. 
We find a stable state for both qqcc and qqbb configurations. Possible decay 
modes of these structures are analyzed. 
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I. INTRODUCTION 



The potentiality of the quark model for hadron physics in the low-energy regime was 
manifest when it was used to classify all the known hadron states. Describing hadrons as 
qq or qqq configurations, their quantum numbers were correctly explained. This assignment 
was based on the comment by Gell-Mann [1] introducing the notion of quark: 'It is assuming 
that the lowest baryon configuration (qqq) gives just the representations 1, 8 and 10, that 
have been observed, while the lowest meson configuration (qq) similarly gives just 1 and 8\ 
Since then, it is assumed that these are the two configurations involved in the description 
of physical hadrons. However, color confinement is also compatible with other multiquark 
structures like the tetraquark qqqq introduced by Jaffe [2]. During the last two decades of the 
past century there appear a number of experimental data that are hardly accommodated in 
the traditional scheme defined by Gell-Mann. Besides, new experiments are being designed 
looking for these unusual structures. 

The possible existence of tetraquarks has been suggested in two different scenarios. The 
first one is a system composed of two light quarks (antiquarks) and two heavy antiquarks 
(quarks). Although such an object is experimentally difficult to produce and also to detect 
[3] it has been argued that for sufficiently large heavy quark mass the tetraquark system 
should be bound [4]. The second scenario corresponds to the scalar mesons, J PC = 0++, 
where a huge amount of experimental data is available nowadays. However, in spite of that 
the situation is not yet conclusive. From the theoretical point of view the scalar sector 
is very particular, because to obtain a positive parity state from a qq pair one needs at 
least one unit of angular momentum. Apparently this costs an energy around 1 GeV, since 
similar meson states (1 ++ and 2 ++ ) lie above 1.2 GeV. However, a q 2 q 2 state can couple to 
J PC = ++ without orbital excitation and, as a consequence, it could be a serious candidate 
to describe the lightest scalar mesons. In this paper we focus our attention on the first 
subject, the heavy-light tetraquarks, whose existence has been recently examined [5] from 
experimental data of the SELEX collaboration [6] . 

The stability of the qqQQ system relies on the mass of the heavy quark [4] . The heavier 
the quark the stronger the short-range Coulomb attraction, in such a way that it could 
play a decisive role to bind the tetraquark system. The binding energy of a 1/r potential is 
proportional to the mass and it must lead to a bound state in the limit of infinite quark mass. 
Moreover the QQ pair brings a small kinetic energy into the system and then contributes 
to stabilize it. On the other hand, heavy-light tetraquarks are ideal systems to study the 
interplay between the different quark interactions because chiral symmetry is spontaneously 
broken in the light sector but it is explicitly broken in the heavy one. 

Heavy-light tetraquarks have been studied in the past in different ways. Carlson et al. [7] 
used a potential derived from the MIT bag model in the Born-Oppenheimer approximation. 
The calculations were done by means of the Green's function Monte Carlo Method. They 
found a J p = 1 + isoscalar bbud bound state and they also concluded that the ccud is not 
bound. A different approach was followed by Manohar and Wise [8] who studied systems 
of two heavy-light (Qq) mesons interacting by a potential determined at long distances 
by a one-pion exchange computed using Chiral Perturbation Theory. They found that for 
Q = b this long range potential may be sufficiently attractive to produce a weakly bound 
two-meson state, although states where Q = c are not bound. In the framework of the 
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nonrelativistic quark potential models, Silvestre-Brac and Semay [9] have studied possible 
tetraquark structures using different parametrizations of the Bhaduri potential [10]. They 
found several bound state candidates in the bbqq and bcqq configurations but not for the 
ccqq structure. A different conclusion is obtained by Pepin et al. [11] using a pseudoscalar 
meson-exchange interaction coming from the breaking of chiral symmetry instead of the 
chromomagnetic potential. Their results indicate that such interaction binds the heavy 
tetraquark systems both for Q = b and Q = c. Therefore, the theoretical situation seems to 
be uncertain depending on whether chromomagnetic or chiral interactions are used. Even 
in the case of the same kind of interaction (one-pion exchange) the results of Refs. [8] and 
[11] are different. 

Such a model dependence of the possible existence of tetraquarks claims for calculations 
with interactions constrained in other sectors. For this purpose, the chiral constituent quark 
model of Ref. [12] is an ideally suited starting point. It is based on the idea that between 
the chiral symmetry breaking scale and the confinement scale, QCD may be formulated for 
the light quark sector as an effective theory of constituent quarks interacting through gluons 
and Goldstone modes associated to the spontaneous breaking of SU (2) chiral symmetry. Its 
parameters have been determined in the description of non-strange two- and three-baryon 
systems and the hadron spectra [12-14]. For the present study one needs to include strange 
and heavy flavors and therefore the model has to be generalized, fixing the new parameters 
in the meson spectra. This interaction will be used to solve the Schrodinger equation for 
the tetraquark system using a variational method. For the spatial wave function we assume 
a linear combination of gaussians and we will consider the two possible color configurations: 
{33} and {66}. 

The paper is organized as follows. In Sec. II we describe the chiral constituent quark 
model. The method and details of the calculation are shown in Sec. III. Section IV is 
devoted to the analysis and discussion of the results. Finally, a summary is presented in 
Sec. V. 

II. THE CHIRAL CONSTITUENT QUARK MODEL 

Since the origin of the quark model hadrons have been considered to be built by con- 
stituent (massive) quarks. Nowadays it is widely recognized that the constituent quark mass, 
very different from the current quark mass of the QCD lagrangian, appears because of the 
spontaneous breaking of the original SU (3)l ® SU(3)r chiral symmetry at some momentum 
scale. The picture of the QCD vacuum as a dilute medium of instantons [15] explains nicely 
such a symmetry breaking, which is the most important nonperturbative phenomenon for 
hadron structure at low energies. Quarks interact with fermionic zero modes of the individ- 
ual instantons in the medium and therefore the propagator of a light quark gets modified 
and quarks acquire a momentum dependent mass which drops to zero for momenta higher 
than the inverse of the average instanton size p. The momentum dependent quark mass acts 
as a natural cutoff of the theory. In the domain of momenta k < 1/p, a simple lagrangian 
invariant under the chiral transformation can be derived as [15] 

L = ^(ry^ - MU l5 )i) (1) 
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where U 15 = exp(i7r a A a 7s// 7r ). n a denotes the pseudoscalar fields (n,Ki,rjs) with i=l,...,4, 
and M is the constituent quark mass. An expression of the constituent quark mass can be 
obtained from the theory, but it also can be parametrized as M(g 2 ) = m q F(q 2 ) with 



F(g 2 ) 



A 2 



A 2 + g 2 



(2) 



where A determines the scale at which chiral symmetry is broken. Even if one does not 
believe in instantons as the microscopic mechanism of spontaneous chiral symmetry breaking, 
one has to admit that once a constituent quark mass is generated by some mechanism 
such quarks inevitably have to interact through Goldstone modes. Whereas the lagrangian 
ipii^^d^ — M)ip is not invariant under chiral rotations, the lagrangian of Eq. (1) is invariant 
since the rotation of the quark fields can be compensated renaming the bosons fields. U 15 
can be expanded in terms of boson fields as, 



1 i ~,5 \ a ^ a 
1 + —7 A 7T 



1 



IX 7T + 



(3) 



The first term generates the quark constituent mass and the second one gives rise to a one- 
boson exchange interaction between quarks. The main contribution of the third term comes 
from the two-pion exchange which will be simulated by means of the one-sigma exchange 
potential. Based on the non-relativistic approximation of the above lagrangian, and making 
use of the physical rj instead the octect one (this is the reason why a mixing angle 9 P appears), 
one can write the following potentials between quarks, 



Vps{n 3 ) = 



9 2 ch 



m 



PS 



A 2 

^PS 



47r Ylmirrij Ap s 



m 



-m PS 



PS 



Y(m PS rn) 



A 3 

^PS 



rn 



PS 



Y(Ap 8rij ) 



(* • *s) E (\ a • AJ) , 



a=l 



(4) 



VsiTij) 



9ch 



A 2 



4tt A 2 



-m s 



Y{m s r %j )-^Y{K s n ) 
ms 



(5) 



9ch 



m 



K 



A 2 



As 



(* • E (\ ? • A") , (6) 



a=4 



V(r-) = & m ^ ^ m 



A 3 

Y(m vrij ) - -^Y(A vrij ) 



mi 



{Bi ■ ffj) cos9 p (\i ■ A|) — sinQp 

(7) 



where g c h = m q /f w , the X's are the SU(3) flavor Gell-Mann matrices, and Y(x) is the 
standard Yukawa function. The chiral coupling constant g c h is related to the nNN coupling 
constant by 
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9ch = /3\ 2 9ttNN m l,d / g x 
4tt \5J 4tt m 2 N { ' 

Proceeding on this way one assumes that flavor SU(3) is an exact symmetry, only broken 
by the different mass of the strange quark, rrii is the quark mass, mps, mx and m v are the 
masses of the SU(3) Goldstone bosons, taken to be their experimental values, and mj is 
taken from the PCAC relation m| ~ m ps + ^ m ud [l^l- m t ne heavy quark sector, chiral 
symmetry is explicitly broken and therefore these interactions will not appear. 

For higher momentum transfer quarks still interact through gluon exchanges. De Rujula 
et al. [17] proposed that gluon exchange between constituent quarks can be described as an 
effective interaction according to the lagrangian 

Lgqq = iV^a s ^G»\ c iP (9) 

where A c are the SU(3) color matrices and is the gluon field. Using a nonrelativistic 
reduction one obtains coulomb and contact potentials that can be parametrized as follows 
[18]: 

VooEinj) = \a s X c i ■ X c j | i- - ^L-Bt ■ a 3 5{r t3 ) } . (10) 

In order to obtain a unified description of light, strange and heavy mesons, a running 
strong coupling constant has to be used [18]. The standard expression for a s (Q 2 ) diverges 
when Q — > Aqcd and therefore the coupling constant has to be frozen at low energies 
[19]. We parametrize this behavior by means of an effective scale dependent strong coupling 
constant similar to the one used in Refs. [20,21] 

where /i is the reduced mass of the qq system and cto, /io and Ao are fitted parameters. 
This equation gives rise to a s ~ 0.54 for the light quark sector, a value consistent with the 
one used in the study of the nonstrange hadron phenomenology [12-14], and it also has an 
appropriate high Q 2 behavior, a s ~ 0.127 at the Z mass [22]. In order to avoid an unbound 
spectrum from below the delta function has to be regularized. Taken into account that for 
a coulombic system the typical size scales with the reduced mass, we use a flavor-dependent 
regularization r (/i) = r /fi, 

°™ - isjoo — • (12) 

The other nonperturbative property of QCD, which cannot be explained by the instanton 
liquid model, is confinement. Up to now it still remains a problem to derive this property 
from QCD in an analytic manner. The only indication we have on the nature of confinement 
is through lattice QCD studies. These calculations show that qq systems are well reproduced 
at short distances by a linear potential. This potential can be physically interpreted in a 
picture in which the quark and the antiquark are linked with a one-dimensional color flux 
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tube or string with a string tension a q g and hence the qq potential is proportional to the 
distance between the quark and the antiquark. However, pair creation screens the potential 
at large distances [23] . A screened potential simulating the results of lattice calculations is 
given by 

Vcon^j) = {~a c (1 - e~^) + A}(A^ • X c j) (13) 

where A is a global constant to fit the origin of energies. At short distances this potential 
presents a linear behavior with an effective confinement strength a = a c fi c (X c i ■ X c j) while 
it becomes constant at large distances. 

Let us resume the different pieces of the interacting potential as a function of the quarks 
involved: 

' (ij) = (qq) Vcon + Voge + V PS + V S + V V 

(ij) = ( qs ) Vcon + V OGE + V s + V k + V v 
Vij = I (ij) = (ss) Vcon + Voce + V s + V v (14) 

(ij) = (qQ) =>- Vcon + Voce 
. (ij) = (QQ) Vcon + Voge 

The corresponding qq potential is obtained from the qq one as detailed in Ref. [24]. In 
the case of V^(fjj), where G-parity is not well defined, the transformation is given by 
\\-\% — > A" • (A2) 71 , which recovers the standard change of sign in the case of the pseudoscalar 
exchange between two nonstrange quarks. Assuming a nonrelativistic expression for the 
kinetic energy, the hamiltonian of the system takes the form: 




tf = E(^ + £r)+£^ ( 15 ) 

where is given by Eq. (14). 



III. TETRAQUARK WAVE FUNCTION 

For the description of the qqQQ system we introduce the Jacobi coordinates 
x = fi-f 2 y = r 3 -f 4 

z = 

mi + m 2 m 3 + m 4 

where particles 1 and 2 are quarks and 3 and 4 antiquarks. The four-body problem is solved 
using a wave function that includes all the possible flavor-spin-color channels that contribute 
to a given configuration. For each channel s, such a wave function will be a tensor product 
of a color (g s ), flavor (f s ), spin (x s ) and spatial (R s ) wave functions 

I <p s >= «? s (1234)/ s (1234) Xs (1234)i? s (1234) (18) 

Concerning the spatial wave function, the most general one with L = may depend on 
six scalar quantities 



R 



Em 



(16) 
(17) 
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.R s (1234) = R s (x ,y ,z ,£ • y,x ■ z,y • z) 



(19) 



Most part of tetraquarks studies have been done by means of a completely symmetric spatial 
wave function, depending only on the quadratic terms, x 2 ,y 2 and z 2 . This is a restricted 
choice which we will analyze in our calculation. We define our variational spatial wave 
function as a linear combination of gaussians 



i? s (1234) = ^^ J R« = E^ ) e- a 



(20) 



i=i 



where n is the number of gaussians we use to expand the spatial wave function of each 
color-spin-flavor component. 

With respect to the color wave function, one can couple the two quarks (1,2) and the 
two antiquarks (3, 4) to color singlet in different ways 



113) 124 > , 

114) I23 > , 
3l2, 334 > , 



813, 824 > 

814, 823 > 

612,634 > 



(21) 
(22) 
(23) 



The first two couplings are convenient for asymptotic meson-meson channels (or meson- 
meson molecules) while the third one is more appropriate for tetraquark bound states. 
With our choice of the Jacobi coordinates the last color basis results to be more suitable 
and essentially to treat the Pauli principle in an easy way. 
The spin part of the wave function can be written as 



S, >= [(12) 5l2 (34) S34 ] 5 



(24) 



where the spin of the two quarks is coupled to and that of the antiquarks to 534. Then 
we will have the following basis vectors as a function of the total spin S=0,l,2: 



S = 




[(12)o(34) ] 
[(12)1(34)4 



(25) 



5 3 >= [(12)0(34)4 

5 4 >= [(12)^34)0]! 

5 5 >= [(12) 1 (34) 1 ] 1 



(26) 



S = 2 



{\S 6 >= [(12)x(34) 



1J2 



(27) 



Concerning the flavor part, we consider light quarks those with flavor u and d and heavy 
ones those with flavor s, c and b. The heavy quarks have isospin zero so they do not 
contribute to the total isospin. Therefore we can classify the tetraquark wave function by 
the isospin of the light quarks I = 0, 1. Taken into account all degrees of freedom, the Pauli 
principle must be satisfied for each subsystem of identical quarks (antiquarks). It imposes 
restrictions on the quantum numbers of the basis states, which is the justification to use the 
coupling [(qq)(QQ)}. 
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Using the wave function described above, we search for a variational solution of the 
hamiltonian of Eq. (15). The spin, color and flavor parts are integrated out and the (3^ 
coefficients of the spatial wave function are obtained by solving the system of linear equations 

n 

E E #° [(RP I H - E (R { ?\Rf)*sA = o v j, s ' (28) 

s i=l 

once the eigenvalues E{a^\ /W) are obtained using a minimization procedure. In practice, 
we use the MINUIT package [25] to obtain the energies. The stable tetraquark states 
are identified by comparing the obtained eigenvalues with the corresponding meson-meson 
threshold calculated with the same hamiltonian. 

IV. RESULTS 

As we have already discussed, the parameters appearing in our hamiltonian correspond- 
ing to the light sector are fixed from the NN interaction [12]. However there are others 
which cannot be determined in this way either because the NN interaction does not depend 
on them (e.g. the confinement parameters) or because they appear as a consequence of the 
generalization to strange and heavy flavors. Among the thoughtful criteria for fixing these 
remaining parameters the most commonly used one is to fit the meson spectra. Being the 
meson masses the ones that will govern the tetraquark thresholds and therefore will deter- 
mine if the tetraquark is bound or not, we think that a good fit of meson spectra must be 
the most important criterium. This criterium is the one proposed long ago by Isgur and 
collaborators [26] and also by Manohar and Wise [8], and it has been also adopted in many 
other works, which use the Bhaduri potential whose parameters are fitted to charmonium 
[9]. Therefore, we have performed fits of the qq sector using the tetraquark wave functions 
in the limit where the two component mesons are isolated, what is equivalent to solve the 
qq system [26]. The complete set of parameters is shown in Table I and the corresponding 
meson ground state masses are given in Table II. 

Tetraquarks will be stable under the strong interaction if their total energy lies below all 
the possible, and allowed, two-meson thresholds. Therefore we will use the quantity 

AE = E T (qqQQ) - E mi (qQ) - E m2 (qQ) (29) 

to discriminate the stable tetraquark states. 

We focus our attention on the lowest L = 0, positive parity, states of the qqQQ system 
with Q = b,c,s and q = u,d for all possible spin-isospin combinations. As a first step we con- 
sider a spatial wave function with a single gaussian and we compare the results obtained for 
the tetraquark energies using the most general wave function [see Eq. (20)] or a completely 
symmetric spatial wave function. The energy difference obtained are always smaller than 
1%. Therefore we will use for the spatial wave function combinations of the quadratic terms, 
R(x 2 , y 2 ,z 2 ), except for those spin-isospin channels which are Pauli forbidden if we use such 
a symmetric spatial wave function. For these channels, which correspond to (S, I) = (0, 0) 
and (2, 0), we make the calculation with the full wave function. 
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The variational calculation is done using as spatial wave function a sum of gaussians. 
We start with one gaussian and we increase its number until convergence in the tetraquark 
energy is reached. The typical number of gaussians needed is six. 

The results are presented in Table III. Our calculation predicts two bound states, one for 
the qqbb system and one for the qqcc, in both cases with (S,I) = (1,0). No bound states are 
found for the qqss system, the (S,I) = (1,0) channel being again the most attractive one. 
The other (S, I) channels are strongly repulsive, specially those where the Pauli principle 
forbids the spatially symmetric wave function, except for the (S, I) = (1, 1) and (2, 1) for 
the qqbb system. Comparing with the existing literature we find two different types of 
calculations: those including only confinement and one-gluon exchange, Refs. [9] and [28] 
and those considering different models of Goldstone boson exchanges between quarks, Refs. 
[8] and [11]. Models including only one-gluon exchange predict systematically a bound state 
in the (S,I) = (1,0) channel for the qqbb system and no bound state for the qqcc system. 
Among the others one finds different conclusions. In Ref. [8] a weakly bound state is found 
only in the qqbb system, whereas Ref. [11] reports bound states for both qqbb and qqcc. 

To clarify the importance of the different pieces of the interacting hamiltonian we have 
redone our calculations switching-off the pseudoscalar part, Eq.(4), of the Goldstone boson 
exchanges. The results are shown in Table IV. A first glance to the results tells us that the 
bound states have disappeared. This part of the interaction plays a relevant role to bind 
the (S,I) = (1,0) channel [11] and it also favors the binding in the (S,I) = (0,0) but not 
in the other channels. One can easily understand the difference with the results of Ref. [9] , 
quoted in Table IV, due to the larger strong coupling constant used in these works. This is 
due to an obvious renormalization of the gluon exchange parameters in order to reproduce 
the meson spectra in two different formalisms (gluon exchange alone and gluon exchange + 
boson exchanges). 

There is an important aspect that should be emphasized at this point. For the predicted 
bound systems the thresholds are determined by D and B mesons. Such mesons are described 
just in terms of a confining and a one-gluon exchange interaction [see Eq. (14)]. Therefore, 
one could use the same restricted interaction to calculate the total energy of the tetraquark, 
but one would obtain different results (see Table IV). The reason for this discrepancy stems 
from the fact that in this case there also appear interactions between light quarks that need 
to be described by means of a more elaborated interaction. As we have previously mentioned 
this shows the importance of reproducing the full meson spectrum (and not a reduced set 
of states) to make thoughtful predictions. 

The structure of our interaction also allows for a study of the influence of the different 
color configurations in the tetraquarks binding energy. In Table V we present the results 
for the probabilities of the two color components {33} and {66} in the tetraquark wave 
function. We observe that for the bound states channel the probability of {66} is almost 
negligible, but its influence increases when the heavy quark mass decreases. Similar results 
have also been reported in the context of a pure one-gluon exchange model [4]. However, 
for unbound channels this probability tends to increase reaching values of 25% for the qqss 
(S, I) = (0, 1). This would be an indication that the {66} channel may be important for the 
pure light tetraquark sector and therefore should not be neglected a priori. 

The dominance of the {33} channel has a more clear physical interpretation if we change 
the color basis from {| 3i 2 , 3 34 >, | 6i 2 , 6 34 >} to {| li 3 , 1 24 >, | 813, 8 24 >}. As stated above, 
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the second basis is more appropriate to describe asymptotic channels because the first vector 
describes two physical mesons, whereas the second one describes two colored meson states. 
The relation between them is given by 

I 3i2, 334 > 

I 612, 634 > 

Being the | 3i2, 334 > channel the dominant one, the contribution to the tetraquark binding 
energy of the colored meson channel is more important than the two physical meson states. 
This fact can be interpreted in the sense that our tetraquarks are a more complicated object 
than a pure meson-meson resonance. 

Let us finally discuss the possible decay channels of these structures. The allowed decay 
modes depend on the relationship between the tetraquark mass and the sum of the masses 
of the possible decay products. In the case of the charmed tetraquark the results shown in 
Table III satisfy mp + mp* > m( qq )(cc) > 2mc, and therefore the strong decay is forbidden, 
being the most probable decay the electromagnetic process T[(qq)(cc)] — > D + D°'y, that 
could be identified detecting a relatively soft photon, E 1 < 140 MeV. In the case of the 
qqbb structure we obtain (see Table III) TTi(q q )(bb) ~~ 2m# < 0, and therefore it can only decay 
through two consecutive weak processes. 

V. SUMMARY 

We were faced in this work with the possible existence of tetraquark structures with 
heavy-light flavors in a variational approach. To be as much predictive as possible we have 
started from a chiral constituent quark model which is able to describe the NN interaction 
and also gives a reasonable description of the non-strange baryon spectrum. The gener- 
alization to strange and heavy flavors has been asked to produce a nice fit of the meson 
spectrum. We have proved that for the tetraquark description it is enough to consider a 
simplified variational wave function depending only on the square of the Jacobi coordinates 
except for those cases for which this wave function is forbidden by the Pauli principle. We 
have found that the influence of the {66} color channel could be important in the light sec- 
tor. For the heavy flavors its contribution is almost negligible which suggests a non-trivial 
color structure for the tetraquark in terms of meson-meson components. 

Our calculation predicts the existence of only one bound tetraquark with (S, I) = (1, 0) 
in both the bottom and the charm sectors. While the first one should decay weakly, the 
second one would do it electromagnetically. 
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TABLES 



TABLE I. Quark-model parameters. 



a c (MeV) 


430.0 


ms (fm x ) 


3.42 


m u4 (MeV) 


313 


He (fm" 1 ) 


0.7 


m PS (far 1 ) 


0.7 


m s (MeV) 


555 


9lJ^ 


0.54 


rriK (fm -1 ) 


2.509 


m c (MeV) 


1752 


a 


2.1181 


m v (fm -1 ) 


2.772 


m b (MeV) 


5100 


A (fm- 1 ) 


0.113178 


A P5i5 (fm- 1 ) 


4.2 


Op (rad) 


-0.2618 


Mo (MeV) 


36.9762 


A KtV (fm" 1 ) 


5.2 


f (fmMeV) 


28.17 






A (MeV) 


181.1 







TABLE II. Meson spectra for the model described on the text. Experimental data are taken 
from Ref. [27]. All masses are given in MeV. 



Meson 


7T 


V 


P 


rf 


4> 


K 


K* 


Result 


139.5 


571.6 


771.6 


955.7 


1020.4 


496.0 


910.5 


Experiment 


138.04 


547.30 


771.10 


957.78 


1019.46 


495.0 


891.66 


Meson 


D 


D* 


Vc 


JM> 


B 


B* 


T 


Result 


1883.2 


2008.9 


2989.9 


3097.0 


5280.8 


5321.2 


9505.0 


Experiment 


1867.7 


2008.9 


2979.7 


3096.87 


5279.2 


5325.0 


9460.30 



TABLE III. E T (qqQQ) and AE energies in MeV. 







(0,0) 


(0,1) 


(1,0) 


(1,1) 


(2,0) 


(2T) 


[(«)(«)] 


Ex 


2239 


1804 


1528 


1940 


2681 


2020 




AE 


+ 1247 


+812 


+126 


+538 


+860 


+199 


[(qq)(cc)} 


Ej> 


4351 


4150 


3764 


4186 


4849 


4211 




AE 


+585 


+384 


-129 


+293 


+830 


+192 


[(qq)(bb)} 


Ex 


10820 


10690 


10261 


10698 


11350 


10707 




AE 


+258 


+128 


-341 


+96 


+708 


+65 
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TABLE IV. Ex(qqQQ) and AE in MeV neglecting the pseudoscalar part of the Goldstone 
boson exchanges. 



(S,I) 




(0,0) 


(0,1) 


(1,0) 


(1,1) 


(2,0) 


(2,1) 


[(qq)(ss)\ 




2630 


1999 


1996 


2139 


2772 


2186 




AE 


+1638 


+1007 


+594 


+737 


+951 


+365 


[(qq)(cc)] 


E T 


4666 


4175 


4101 


4231 


4851 


4254 




AE 


+900 


+409 


+208 


+338 


+832 


+235 




AE from Ref. [9] 






+19 








[(qq)(bb)} 


Ef 


11136 


10735 


10612 


10743 


11352 


10752 




AE 


+574 


+173 


+10 


+141 


+710 


+110 




AE from Ref. [9] 






-131 


+56 




+30 



TABLE V. Probability of the different color channels. 



(s,i) 


(0,1) 




(1,0) 






P {33} 


P {66} 


P {33} 


P {66} 


[(qq)(ss)} 


0.7469 


0.2531 


0.9706 


0.0294 


[(qq)(cc)} 


0.8351 


0.1649 


0.9940 


0.060 


Km) (55)] 


0.9895 


0.0105 


0.9993 


0.0007 
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